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Abstract 

!>■ ' 

pq ' The production and propagation of mesons (it, r], p,u,$>, K, K, J/^f) in proton- 

nucleus and nucleus-nucleus collisions from 1 - 200 GeV/u is studied within the 
covariant transport approach HSD, which explicitly allows to investigate selfen- 
ergy effects of the hadrons at finite baryon density. Whereas the experimental 
pion and K + spectra can be described without introducing any selfenergies for 
the mesons, the K~ yield in Ni + Ni collisions is underestimated by a factor of 

^^ \ 5-7 at 1.66 and 1.85 GeV/u. However, introducing density dependent antikaon 

masses in line with effective chiral Lagrangians a satisfactory agreement with 
the data is achieved. A dropping of the p-meson mass with baryon density, as 
suggested by QCD sumrule studies, is proposed to explain the dilepton spectra 
for S + Au and Pb + Au at SPS energies, which indicates independently that 
a partial restoration of chiral symmetry might be found already in the present 



OS 



o 

experiments. 



1 Introduction 



Ctf 



The study of hot and dense nuclear matter via relativistic nucleus-nucleus collisions is 
the major aim of high energy heavy-ion physics. Nowadays, the search for a restoration 
of chiral symmetry at high baryon density or for a phase transition to the quark-gluon 
plasma (QGP) is of specific interest. However, any conclusions about the hadron 
properties at high temperature or baryon densities must rely on the comparison of 
experimental data with theoretical approaches based on nonequilibrium kinetic theory. 
Among these, the covariant RBUU approach [|TJ, g, the QMD |§ or RQMD model |§, 
and the HSD approach || have been successfully used in the past. As a genuine feature 
of transport theories there are two essential ingredients: i.e. the baryon (and meson) 
scalar and vector selfenergies - which are neglected in many approaches - as well as 
in-medium elastic and inelastic cross sections for all hadrons involved. 

Selfenergy effects in the production of particles have been found so far for antipro- 
tons at SIS energies [|, [|J, though the actual magnitude of the attractive p-potential 
in the nuclear medium is still a matter of debate since the strong p annihilation with 
nucleons is not sufficiently under control. As advocated in Refs. || ||, [U], |ll| on the 



basis of effective chiral Lagrangians also antikaons should feel strong attractive forces 
in the medium so that their production threshold should be reduced at finite baryon 
density. The vector mesons p and u, furthermore, are expected to drop in mass with 



baryon density according to QCD sumrule studies [jT2j , which is basically a consequence 
of the dropping scalar quark condensate < qq > with the quark density < q^q > |Tj| . 
Whereas a direct enhancement of the p-meson yield is hard to observe experimentally 
due to the short lifetime of the p and the strong final state interactions of the pions, 
dilepton spectroscopy is expected to provide valuable information on the p spectral 
function in the dense medium. 

The present article is organized as follows: Section 2 contains a brief description 
of the transport approach employed as well as a specification of the meson selfenergies 
incorporated in the calculation. Section 3 is devoted to a presentation of the calculated 
results for ir,K + , and K~ spectra in comparison to available data at SIS, AGS and 
SPS energies. Section 4 concentrates on dilepton physics at SPS energies while Section 
5 concludes this study with a summary and discussion of open problems. 

2 Ingredients of the transport approach 

In this work the dynamical analysis of p+A and A+A reactions is performed within the 
HSD [] approach which is based on a coupled set of covariant transport equations 
for the phase-space distributions fh{x,p) of hadron h ]|j, ||, i.e. 

{(n M - n^u" - M* h diu s h ) % + (n^x + m*^u*) %} f h (x, P ) 
J2 Jd2d3dA.. . [G t G]i2-+34...<J 4 (n + n 2 - n 3 - n 4 . . .) 

ih-ihi... 
X {fh 3 (x,P3)fh 4 ( x >P4)fh(x,P)fh 2 (x,P2) 

- fh(x } p)f h2 (x,p 2 )f h . i (x,p 3 )f h4 (x,p 4 )}... . (1) 

In Eq. ([[]) Uft(x,p) and Uj*(x,p) denote the real part of the scalar and vector hadron 
selfenergies, respectively, while [G + G]i2^34...£ 4 (II + n 2 — n 3 — n 4 . . .) is the 'transition 
rate' for the process 1 + 2— > 3 + 4 + ... which is taken to be on-shell in the semiclassical 
limit adopted. The hadron quasi-particle properties in ([!]) are defined via the mass-shell 
constraint 0, 

5(n M iP-M* 2 ) , (2) 

with effective masses and momenta (in local Thomas-Fermi approximation) given by 

M*(x,p) = M h + U^(x,p) 

U.»(x,p)=if-Ui!(x,p) , (3) 



h2hzh±. 



while the phase-space factors 



f h {x,p) = l±f h {x,p) (4) 



1 Hadron String Dynamics 



are responsible for fermion Pauli-blocking or Bose enhancement, respectively, depend- 
ing on the type of hadron in the final/initial channel. The dots in Eq. (]l|) stand 
for further contributions to the collision term with more than two hadrons in the fi- 
nal/initial channels. The transport approach flU) is fully specified by U^(x,p) and 
U£(x,p) (/i = 0,1,2,3), which determine the mean-field propagation of the hadrons, 
and by the transition rates G^G<5 4 (. . .) in the collision term, that describes the scat- 
tering and hadron production/absorption rates. 

The scalar and vector mean fields U^ and Ufi for baryons are taken from Ref. |J 
and don't have to be specified here again. The pions and 77's as Goldstone bosons are 
expected not to change their properties in the medium; they will be treated as bare 
particles throughout all calculations. Furthermore, the K + meson energy changes only 



very moderately with baryon density according to Kaplan and Nelson [10, 14, [TjJ due 
to a partial cancellation of the scalar and vector kaon selfenergies. Thus, they are also 
produced and propagated as free particles. The antikaon and vector meson potentials 
in the medium, however, have to be specified more explicitly. 

2.1 K~, p,u, $ selfenergies 

As in case of antiprotons there are a couple of models for the antikaon and vector meson 
selfenergies which differ in the actual magnitude of the meson potential. Without going 
into a detailed discussion of the various approaches we adopt the more practical point 
of view, that the actual K~,p,u and $ selfenergies are unknown and as a guide for 
our analysis use a linear extrapolation of the form (for the meson x), 

m t(PB) = ra° I 1 - a x — > m q + m ? , (5) 

V PoJ 

with a x ~ 0.2 for antikaons, which leads to a fairly reasonable reproduction of the 



antikaon mass from Refs. fLOl , [L4| and the recent results from Waas, Kaiser and Weise 



TI[, a x ~ 0.18 for p and u mesons according to Hatsuda and Lee [0 and a x ~ 0.025 
for the $ meson. We note that the dropping of the meson masses is associated with 
a corresponding scalar energy density in the baryon/meson Lagrangian, such that the 
total energy- momentum is conserved during the heavy-ion collision (cf. Ref. ||). 

2.2 Elastic and inelastic reaction channels 

Baryon-baryon (BB) collisions are described using free differential cross sections from 
Ref. [|1| for invariant energies ^/s < 2.6 GeV and by the LUND string formation 



and fragmentation model [0 for ^fs > 2.6 GeV, which generates the hadronic final 
states of a BB collision. The same concept is used for meson-baryon (mB) reactions, 
where for ^Js < 2.3 GeV differential cross sections from Ref. [jnj are employed whereas 
the LUND model is used for higher invariant energies. Meson-meson (mm) reactions 
(e.g. mi — > p, 7rp — > $,7rp — > a\) are described within the Breit-Wigner resonance 



picture using branching ratios from the nuclear data tables fl8 |. In all reaction chan 



nels the thresholds are shifted according to the actual mass of the hadrons (at finite 
baryon density ps = Jjfij^, where j^x) is the local baryon current). Also within the 



Breit-Wigner resonance formation the actual masses of the hadrons are used, whereas 
their width is corrected according to the local phase space for the decay. For the de- 
tailed parametrizations employed (including also baryon-hyperon and meson-hyperon 
channels) the reader is refered to Refs. ]5|, [19], ^0] . 



2.3 Optimizing for high baryon density 



Pb+Pb, central collision 



10 



S 8 



of 6 
■ 

a 
Q. 

r 4 

a 











^•^ r mm 


= 2P 
p 

1 mi 


n= 3 Po 


P = 

r mm 

P, 


^^^^^ 


=5p — — 

nin "0 



10 



100 



E, ab (GeV/u) 



Figure 1: The quantity F ([]) for central Pb + Pb collisions as a function of the 
bombarding energy per nucleon for 4 different cuts in p m i n - 

In order to probe the restoration of chiral symmetry at high baryon density in 
nucleus-nucleus collisions, one has to perform experiments with heavy nuclei (e.g. Pb 
+ Pb) and optimize the beam energy to achieve a large volume of high baryon density 
for a sufficiently long time. In this respect central collisions of Pb + Pb have been 
investigated within the transport approach specified above and the 'stopped' baryon 
density pl(t) - including only baryons with rapidity \y\ < 0.7 in the cms - has been 
computed in a central volume V = ttR^/'Jcm with R = 4 fm, while Jcm is the Lorentz 
factor in the nucleus-nucleus center-of-mass system. Since we are interested in high 
baryon densities above some value p m i n for long times, we consider the quantity 



F = J dt (p s B (t) - Pmin ) Q(p s B (t) - p„ 



(6) 



which should serve as a useful guide in the optimization problem. The quantity F (|6]) 
is displayed in Fig. 1 for central collisions of Pb + Pb from 1 - 200 GeV/u for different 
values of p min from 2 - 5 p (Po ~0.168/m -3 ). Thus optimal bombarding energies for 
baryon densities above 4p should be around 20 - 30 GeV/u in order to explore the 
properties of an intermediate phase, where the chiral symmetry might approximately 
be restored and the hadron masses (except for the Goldstone bosons) might be close 
to their current quark masses m q + irtq. 



3 Meson production in nucleus-nucleus collisions 



3.1 Pions 

The pions as the lightest Goldstone bosons are not expected to change their properties 
in the dense nuclear medium significantly - except in a quark-gluon-plasma (QGP) 
phase - such that their production and propagation should be reasonably described 
without introducing any selfenergies. As an example for SIS energies we show in Fig. 2 
the calculated results || for transverse it spectra in Ar + Ca collisions at 1.5 GeV/u 
(for 0.68 < yiab < 0.84) in comparison to the data of the TAPS collaboration |21| (open 
squares), which are described over four orders of magnitude with relative deviations of 
less than 30%. Similar results have been obtained for various systems at 1 - 2 GeV/u 
in Ref. ||19||. Quantitatively similar experiences have been made at AGS energies as can 
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Figure 2: The calculated transverse ir° spectra (solid line) in Ar + Ca collisions at 1.5 
GeV/u (for 0.68 < y\ a b < 0.84) in comparison to the data of the TAPS collaboration [pl| 
(open squares) 
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Figure 3: The calculated transverse mass spectra for tt + (solid lines) in ntral collisions 
of Si + Al at 14.6 GeV/u for different rapidities in laboratory (y = 0.9,1.7,2.7) in 
comparison to the data from Ref. 



be extracted from Fig. 3 where the transverse n + mass spectra for central collisions of 
Si + Al at 14.6 GeV/u are displayed for 3 different rapidities in the laboratory system 
in comparison to the data from Ref. |22|j . 

As an example for the pion yield at SPS energies we show in Fig. 4 the ir~ rapidity 
distribution for central S + S collisions at 200 GeV/u in comparison to the data from 
Ref. |23| (open squares), which are approximately of Gaussian shape. At midrapidity 
(y = 0) here the 7r~ density is about a factor of 7 higher than the proton density. This 
indicates that the available energy is dominantly used for mass production (in form 
of pions) and that during the longitudinal expansion of the 'hadronic fireball' meson- 
meson reaction channels should occur more frequent than meson-baryon or baryon- 
baryon reactions. 
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Figure 4: The n~ rapidity distribution (solid line) for central S + S collisions at 200 
GeV/u in comparison to the data from Ref. |23j (open squares). The dotted line shows 
the calculated proton rapidity distribution in comparison to the respective data from 
Ref. p3f (full squares). 
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Figure 5: The preliminary rapidity distribution of negative hadrons from NA49 [M 
(full squares) in comparison to the HSD results p6 |. The solid line corresponds to a 
calculation (at b = 2fm) including the dropping meson masses from Eq. (|5]), whereas 
the dotted line results from a calculation with bare meson masses. The open squares 
are obtained by reflecting the full squares at midrapidity. 
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The shape of the pion rapidity distribution is not changed significantly when going 
over to central collisions of Pb + Pb at 160 GeV/u as can seen from Fig. 5, where 
the preliminary rapidity distribution of negative hadrons (essentially n~ , K~ and p) 
from NA49 |J4j is shown in comparison to the HSD results p6| . Here the solid line 



corresponds to a calculation (at b = 2fm) including the dropping meson masses from 
Eq. (HD, whereas the dotted line results from a calculation with bare meson masses. 
The broadening of the rapidity distribution around midrapidity (y « 3) in the dropping 
mass scenario is due to pions from p and to decays, which are produced with a wider 
distribution in rapidity in this case. Since the proton distribution dN p /dy ^40 — 45 at 
midrapidity (cf. Fig. 22 in Ref. ||) the tt~ to proton ratio is only about 4 for central 
Pb + Pb collisions such that the baryon density in the expanding 'hadronic fireball' is 
significantly higher than that for central S + S collisions. 

3.2 Kaons and antikaons 

Whereas kaons should feel a slightly repulsive potential in the nuclear medium ac- 
cording to the approach by Kaplan and Nelson [|K| or Waas, Kaiser and Weise |TT| , 
the antikaons should experience a stronger attractive potential at finite baryon density, 



which is also supported by K atomic data [27 . As a first order approximation we thus 



assume the K + potential or selfenergie to be zero and adopt the linear parametrization 
for the in-medium K~ mass from Eq. (|5]). For the detailed reaction channels and cross 
sections considered the reader is refered to Ref. [^(| . 



The l.h.s. of Fig. 6 shows the calculated results for the inclusive K + invariant 
cross section for Ni + Ni collisions at 0.8, 1.0 and 1.8 GeV/u at 9i a t, =44°, that have 
been transformed to the nucleus-nucleus cms, in comparison to the preliminary K + 



spectra from Ref. |28], |29fl . Since the data can be described quite reasonably at all 
energies from 0.8 - 1.8 GeV/u, apparently no selfenergy effects are needed for K + 
mesons. This finding is also in accordance with earlier studies on K + production in 



nucleus- nucleus fl30| £31], [3J| and proton- nucleus collisions p3fl . On the hand, due to 
the rather stable quasi-particle properties of the kaons in the medium, they qualify as 
probes in connection with the nuclear-equation-of-state (EOS) as suggested early by 
Aichelin and Ko [[34 . 

The r.h.s. of Fig. 6 shows the calculated K~ spectra for Ni + Ni at 1.85 GeV/u 
at 0° with respect to the beam axis in the nucleus-nucleus cms in comparison with 
the data of Ref. |35[ (full squares) and the preliminary data for Ni + Ni at 1.8 GeV/u 



from Ref. [|3(| (open dots). The dashed line reflects a calculation including the bare 
antikaon mass without any antikaon absorption, while the dash-dotted line includes 
antikaon absorption, which reduces the cross section on average by a factor of 5 for 
the Ni + Ni system. However, the data are underestimated severely in the bare K~ 
mass approximation. The solid line in Fig. 6 (r.h.s.) shows the result of a calculation, 
where the K~ mass drops with baryon density according to Eq. (^) with a x = 0.2 
including also antikaon reabsorptionP]. With increasing a x not only the magnitude 

2 For practical purposes one should consider a x to be a free parameter to be fixed in comparison to 
the experimental data in order to learn about the magnitude of the antikaon selfenergy. In fact, we 
obtain a much better reproduction of the spectra for a x « 0.24, but due to the uncertainties involved 
in the elementary BB production cross sections one cannot determine this value very reliably. 



of the spectrum is increased, but also the slope becomes softer. For a x ~ 0.2 we 
still underestimate the experimental spectra slightly, but it is clearly seen that quite 
sizeable antikaon attractive selfenergies are needed to reproduce the data. This finding 



is also in line with an independent calculation by Li, Ko and Fang |37 
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Figure 6: (l.h.s.) The calculated results for the inclusive K + invariant cross section for 
Ni + Ni collisions at 0.8, 1.0 and 1.8 GeV/u at 6i a b =44°, that have been transformed 
to the nucleus-nucleus cms, in comparison to the preliminary K + spectra from Ref. [^, 
|2"9fl . (r.h.s.) The calculated K~ spectra for Ni + Ni at 1.85 GeV/u at 0° with respect 
to the beam axis in the nucleus-nucleus cms in comparison with the data of Ref. [^ 
(full squares) and the preliminary data for Ni + Ni at 1.8 GeV/u from Ref. [|36| (open 
dots). The dashed line reflects a calculation including the bare K~ mass without any 
antikaon absorption, while the dash-dotted line includes antikaon absorption. The solid 
line corresponds to a calculation with a dropping antikaon mass according to Eq. (|5]) 
for a x = 0.2. 

The enhanced production of strangeness is also known from experiments at AGS 
and SPS energies |3£|. As an example Fig. 7 shows the measured K + /n + ratio for pp, 
Si + Al, Si + Cu, Si + Au and Au + Au at AGS energies 
a factor of 3 with the number of participating nucleons. 
with in-medium meson masses (solid line) 
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which increases by about 
Whereas a HSD calculation 
approximately reproduces this trend, 
the same calculation with bare meson masses (lower dotted line) underestimates the 
K + /tt + ratio significantly. The actual enhancement in the dropping mass scenario is 
due to a large contribution from the channel meson + meson — > K + K~ , which is 
enhanced sizeably for a dropping antikaon mass. Differential transverse momentum 
spectra for kaons and antikaons at midrapidity for central Au + Au collisions at 10.8 
GeV/u should shed some further light on this issue. Note, that RQMD calculations 
for these systems also underestimate the K + /ix + ratio by about a factor of 1.5. 
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Figure 7: The measured K + /n + ratio for pp, Si + Al, Si + Cu, Si + Au and Au + 
Au at AGS energies |39[ (full dots) in comparison to HSD calculations with in-medium 
meson masses (solid line) || and bare meson masses (lower dotted line). 

4 Electromagnetic probes 

Photons and dileptons are particularly well suited for an investigation of the violent 
phases of a high-energy heavy-ion collision because they can leave the reaction volume 
essentially undistorted by final-state interactions. Whereas the signal from direct pho- 
tons is largely covered by the electromagnetic decays of light neutral mesons (7r°,?7), 
e + e~ or p + p~ pairs at higher invariant masses do not suffer that much from large 
background contributions. Indeed, dileptons from heavy-ion collisions have been ob- 
served by the DLS collaboration at the BEVALAC gl], |42|, 
HELIOS |5], g§, NA38 gTfl and NA50 [g| collaborations at SPS energies 



and by the CERES [44 



4.1 



e T e pairs 



Quite some years ago it has been found within microscopic transport studies at BE- 
VALAC/SIS energies (49| that above about 0.5 GeV of invariant mass (of the lepton 
pair) the dominant production channel is from 7r + 7r~ annihilation, such that the prop- 
erties of the short lived p meson could be explored at high baryon density. The data 
available so far, however, did not allow for a closer distinction of the various models 
proposed. At SPS energies the enhancement of the low mass dimuon yield in S + W 
compared to p + W collisions ]45[ has been first suggested by Koch et al. [|5(J to be due 
to 7r + 7r~ annihilation. Furthermore, Li et a/.[^T[ have proposed that the enhancement 
of the e + e~ yield in S + Au collisions - as observed by the CERES collaboration J44] - 
should be due to an enhanced p-meson production (via 7r + 7r~ annihilation) and a drop- 
ping p-mass in the medium. In fact, their analysis - which was based on an expanding 
fireball scenario in chemical equilibrium - could be confirmed within the HSD trans- 
port calculations in Ref. |5I|. However, also a more conventional approach including 
the increase of the p-meson width in the medium due to the coupling of the p, ir, A 



and nucleoli dynamics |53] was found to be roughly compatible with the CERES data. 
On the other hand, the dimuon data of the HELIOS-3 collaboration [45| could only be 
described satisfactorily when including dropping meson masses |3|]. In the following 
some more recent results from Refs. p5|, |^] are reported, where systematic studies on 
the various dilepton channels from 10 - 200 GeV/u have been performed. 
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Figure 8: The calculated dilepton spectra (full solid line) for p + Be at 450 GeV/u 
in comparison with the data from Ref. f33| . The thin lines indicate the individual 
contributions from the different production channels including the CERES-acceptance 
and mass resolution. 

As an example for dilepton spectra at SPS energies Fig. 8 shows the spectral 
decomposition as a function of the e + e~ invariant mass M for p + Be at 450 GeV/c 
in comparison to the data of the CERES collaboration |l|J. In this case the spectrum 



can be fully accounted for by the electromagnetic decays of the r], rf and vector mesons 
p°,uj and $. Contributions from meson-meson channels (tt + 7t~, K + K~, Tip) are of 
minor importance. 

The situation changes quite dramatically when going over to nucleus-nucleus colli- 
sions. For Pb + Au at 160 GeV/u (and semicentral collisions) the dominant yield for 
invariant masses 0.3 GeV < M < 0.7 GeV stems from 7r + 7r~ annihilation (cf. Fig. 9). 
Also in the $ mass regime about 1 GeV there is a large contribution from K + K~ and 
up annihilation to dileptons for both scenarios: with bare meson masses (upper part 
of Fig. 9) and with in-medium meson masses (lower part of Fig. 9). Whereas most of 
the processes (Dalitz and direct decays) occur in the vacuum at zero baryon density, 
the 7T7T — » p° — » e + e~ and direct p° (from BB and mB collisions) decay still occur at 
finite baryon density such that a dropping p mass also leads to a shift of the respective 
contribution to lower invariant masses M. In Fig. 9 both scenarios are compared to 
the preliminary data of the CERES collaboration 55| including the experimental cuts 
in rapidity y, transverse momentum of the leptons as well as the CERES mass reso- 
lution; due to the present statistics, however, there is no unique conclusion since the 
calculation with bare meson masses (upper part) also describes the data except for one 
point at 0.6 GeV. On the other hand, the present preliminary data match well with the 
calculation including the in-medium meson masses. Apart from better statistics also a 
higher mass resolution of the CERES detector (especially in the p, u, $ mass regime) 
should allow to disentangle the different scenarios in the next years. 
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Figure 9: Dilepton invariant mass spectra for semi central collisions of Pb + Au at 
160 GeV/u (full solid lines) in comparison to the preliminary data of the CERES 
collaboration ||55| . The upper part shows the results of a calculation with bare meson 
masses whereas the lower part includes the dropping meson masses (0). 



5 Summary 

In this article the production of secondary particles in proton-nucleus and nucleus- 
nucleus collisions from 1 - 200 GeV/u has been investigated within the covariant trans- 
port approach HSD [J5J. The analysis shows that ir and K + spectra are reasonably 
well described in this energy regime without introducing any medium modifications 
for these mesons (cf. also Ref. |B| in case of pions). This experience is fully in line 
with earlier studies on this subject and the results from independent groups |56], |57| . 
The antikaon spectra, however, are underestimated severely when incorporating only 
bare kaon masses roughly in line with the study by Li et al. ||37|| . When including an 
attractive antikaon potential comparable to that proposed by Waas, Kaiser and Weise 
Tlfl , a satisfactory description of the K~ spectra can be given, both in the actual 
magnitude as well as in the slope. It is worth noting that the 7r-hyperon — > K~N 
production channels play a sizeable role in case of the vacuum antikaon mass, whereas 
their contribution in the 'dropping mass scenario' becomes of minor importance. 

The 'observed' dropping of the antikaon mass with baryon density may be inter- 
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preted as a step towards a partial restoration of chiral symmetry that can already be 
seen at SIS energies (1-2 GeV/u). Similar observations have been made at AGS ener- 
gies (10 - 15 GeV/u) (cf. Fig. 7) as well as at SPS energies (160 - 450 GeV/u), where 
especially a dropping of the p-mass can be used to accurately describe the dilepton 
spectra from heavy-ion reactions |S^, [H], |5I| at invariant masses 0.3 GeV < M < 0.7 
GeV. The enhancement of dileptons (about a factor of 3) seen in the $ mass region is 
essentially due to the meson-meson production channels irp — > $ and K + K~ — ► $. 

Though there are quite a number of indications for dropping meson masses in 
the medium by now, one has to properly examine the possibility that conventional 
many-body effects such as 'resonance-hole' loops |58 may also account for the spectra 



observed. In addition, detailed experimental studies on K~,p and u> production in 
proton (pion) - nucleus reactions should be performed close to threshold energies since 
a 20% reduction of their mass at density po should clearly be visible in the respective 
spectra 
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